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Abstract. Most biological beams bend and twist rela- 
tively easily compared to human-made structures. This pa- 
per investigates flexibility in 57 diverse biological beams in 
an effort to identify common patterns in the relationship 
between flexural stiffness and torsional stiffness. The pat- 
terns are investigated by mapping both ideal and biological 
beams into a meehanospace defined by flexural and tor- 
sional stiffness. The distribution of biological beams is not 
random, but is generally limited to particular regions of the 
meehanospace. Biological beams that are stiff in bending 
are stiff in torsion, while those that bend easily also twist 
easily. Unoccupied regions of the mechanospace represent 
rare combinations of mechanical properties, without prov- 
ing that they are impossible. The mechanical properties of 
biological beams closely resemble theoretical expectations 
for ideal beams. Both distributions are potentially being 
driven by the interdependence of the material and structural 
properties determining stiffness. The meehanospace can be 
used as a broadly comparative tool to highlight systems that 
fall outside the general pattern observed in this study. These 
outlying beams may be of particular interest to both biolo- 
gists and engineers due to either material or structural 
Innovations. 


Introduction 


Flexibility, or the ability to deform in response to a load, 
is a property of most Liological beams (Vogel, 1984; 
Denny, 1988), yet the biological consequences of flexibility 
vary widely. In motile ors ms, flexibility permits the 
relative movements of structural elements in response to 
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internal forces generated by museular contractions or hy- 
drostatie pressures. The flexibility of a fish backbone influ- 
ences the mechanical behavior of the body during undula- 
tory swimming (McHenry et al., 1995; Long and Nipper, 
1996), while the flexibility of mammalian backbones has 
been implicated in locomotor differences between species 
(Gal, 1993). In many sessile organisms, flexibility allows 
structures to passively adjust their posture relative to the 
forces experienced (Wainwright et al., 1976; Vogel, 1984). 
Leaf petioles (Vogel, 1989; Niklas, 199]) and herbaceous 
plants (Ennos, 1993; Etnier and Vogel, 2000) reduce flow- 
induced drag by bending or twisting in response to wind, 
and similar drag-reducing mechanisms have been found in 
hydroid colonies (Harvell and LaBarbera, 1985) and anem- 
ones (Koehl, 1977a). Other flexible organisms take advan- 
tage of external forces to passively orient their filter-feeding 
structures in response to ever-changing flow (Wainwright 
and Dillon, 1969; Koehl, 1977b; Harvell and LaBarbera, 
1985: Best, 1988). Thus, the ability to deform in response to 
loads is observed in both motile and sessile organisms living 
in either an aquatic or a terrestrial environment, apparently 
independent of phylogenetic affiliations. Such convergence 
may be viewed as a red flag indicating the tremendous 
importance of flexibility in the design of biological organ- 
isms (Lauder, 1982. Vogel, 1998). 

Flexibility is measured in terms of stiffness, where flex- 
ural stiffness (EJ in N +m?) represents the resistance of a 
beam (a structure that is long relative to its width) to 
bending, and torsional stiffness (GJ in N : m`) represents the 
resistance of a beam to twisting. Flexural stiffness and 
torsional stiffness are composite variables that are influ- 
enced both by material and structural properties (Wain- 
wright et al., 1976). Every beam is characterized by a 
combination of flexural stiffness and torsional stiffness, and 
the relationship between these two variables determines 
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how the beam responds to a given load. The ratio of flexural 
stiffness to torsional stiffness, commonly termed the twist- 
to-bend ratio, has been used as a dimensionless (and, thus, 
size invariant) index describing the relationship between 
these two variables (Niklas, 1992; Vogel, 1992, 1995; Et- 
nier and Vogel, 2000). The twist-to-bend ratio indicates the 
relative resistance of a beam to bending versus twisting. 
More intuitively, a higher twist-to-bend ratio indicates a 
structure that twists more readily than it bends, without 
reference to the magnitude of either variable. 

While flexibility is a common property of a phylogeneti- 
cally diverse group of organisms, are there any common 
patterns or trends in the relationship between flexural stiff- 
ness and torsional stiffness in biological beams? This paper 
investigates such pattems with a mechanospace defined by 
values of flexural and torsional stiffness. The mechano- 
space, similar to Raup’s (1966) classic morphospace, is a 
broadly comparative tool used to visualize the relationships 
between mechanical variables in biological beams. The 
concept is based on the premise that the mechanical prop- 
erties of flexural and torsional stiffness are common to all 
biological beams. Three variations of this mechanospace 
will be used to compare the patterns of flexibility seen in a 
large diversity of biological structures. First, material and 
structural properties will be used in combination to predict, 
on the basis of principles of engineering beam theory, the 
theoretical relationships between bending and twisting in 
ideal beams. Second, experimentally measured values of 
flexural stiffness and torsional stiffness for biological struc- 
tures will be examined within the context of the theoretical 
distribution, Finally, the relative contribution of overall size 
to the mechanical properties of biological beams will be 
explored. The results suggest that the distribution of bio- 
logical beams within the mechanospace is not random, due 
to the interdependence of material and structural properties 
determining stiffness. 


Materials and Methods 
Distribution of ideal beams 


The distribution of ideal beams was determined using 
principles from engineering beam theory. Importantly, this 
distribution is limited to structures built of a single, isotropic 
material (i.e., the material properties are not directionally 
variable), with precise specifications for the cross-sectional 
shape of the beam in question (Roark, 1943). Additionally, 
engineering beam theory stipulates that the material is lin- 
early elastic, and that the beam is straight and does not vary 
in size or shape along its length, nor does the beam undergo 
deflections greater than 10% of total length (Roark, 1943). 
More complex solutions are required for beams that un- 
dergo larger deflections (e.g., Morgan and Cannell, 1987; 
Morgan, 1989). 


Material properties. The material properties influencing 
flexural stiffness and torsional stiffness are Young’s mod- 
ulus (EF in N- mM) and the shear modulus (G in N- m73), 
respectively. Young’s modulus and the shear modulus are 
related to one another by Poisson's ratio (v), which is the 
dimensionless ratio of the induced strain, causing lateral 
contraction of the specimen, to the applied strain, causing 
the specimen to elongate (Vincent, 1990). Poisson’s ratios 
can vary from O to 0.50 for naturally occurring isotropic 
materials. Mollusc shell has a Poisson’s ratio of about 0.10, 
while rubber has a ratio closer to 0.50 (Denny, 1988) and 
materials such as cornstalks have moderate ratios around 
0.23 (Prince and Bradway, 1969). Commonly occurring 
metals have Poisson's ratios between 0.25 and 0.30 (Niklas, 
2): 

For isotropic materials, the shear modulus is related to 
Young's modulus (Roark, 1943) by: 


E 
Ci S (1) 


A Ea 
Thus, the Young`s modulus for a typical isotropic material 
will range from 2 to 3 times its shear modulus as Poisson’s 
ratio varies from 0 to 0.50 (Wainwright et al., 1976; Niklas, 
oe). 

Structural properties. The structural variables influencing 
flexural stiffness and torsional stiffness are the second mo- 
ment of area (J in m*) and the polar moment of area (J in 
mt), respectively. These variables reflect the geometry of a 
cross section of a beam and are influenced by size, shape, 
and orientation (Roark, 1943). The relationship between / 
and J depends upon the cross-sectional shape of the struc- 
ture in question (Roark, 1943). Formulas for the calculation 
of I and J for most simple shapes can be found in any basic 
engineering handbook (e.g., Gere and Timoshenko, 1984). / 
and J are both proportional to radius to the fourth power, 
hence radius is a very strong determinant of stiffness 
(Roark, 1943). For example, for a beam with a circular 
cross-sectional area 


TF m 
I= anda (2) 
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the value of /// is 0.50. Small changes in the cross-sectional 
shape of a beam can greatly influence the values of / and J; 
thus, the value //J for noncircular cross sections can be 
much higher (Table 1). Note that there may be several 
values of //J for a beam with an asymmetric cross-sectional 
shape, depending on its orientation with respect to the 
applied load (Table 1). 

Relationship between flexural and torsional stiffness. 
El and GJ are dependent variables, with Poisson's ratio 
linking the material properties (E and G), and geometry 
linking the structural properties (/ and J). Theoretically, the 
relationship between flexural stiffness and torsional stiff- 
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Table 1 
Theoretical relatiorships between flexural and torsional stiffness 
Ratio of El/GJ 

If» = 0, If v = 0.50, 
-Honal shape HJ then E/G = 2 then E/G = 3 

0.50 1.00 1.50 

Filhpse 4:1 major:minor axis 027 0.54 0.81 

Ellipse 4:1 major:minor axis 4.25 8.50 [275 








The ratio of //J was determined from beam-theory equations (Niklas, t992) and represents the structural 
contributions to stiffness in these beams. For the ellipses, / is calculated about the axis indicated by the dashed 
line. The ratio of E/G represents the material contribution to beam stiffness as Poisson's ratio was allowed to 
vary from 0 to 0.50. The chosen cross-sectional shapes are broadly representative of biological beams 


(Wainwright, 1988). 


ness can vary only slightly, based on this interdependence. 
The relationships between flexural stiffness and torsional 
stiffness can be calculated for beams of different shapes and 
materials (Table 1). The ratio of £7/GJ for a given circular 
beam can range between | and 1.5 as Poisson's ratio varies 
from 0 to 0.5. In contrast, an elliptical beam with a major to 
minor axis ratio of 4.1 has a range of El/GJ (with EI 
measured about the minor axis) from 8.50 to 12.75 for the 
same range in Poisson's ratio (Niklas, 1992). These partic- 
ular cross-sectional shapes are presented because they are 
broadly representative of biological beams (Wainwright, 
1988). Note that two beams can have the same E//G/J ratio 
despite vast differences in size. shape, and material because 
the ratio is determined by the relative magnitudes of flexural 
stiffness to torsional stiffness, not by their absolute values. 

The relationship between E/ and GJ for ideal beams can 
be mapped into a mechanospace defined by these variables 
(Fig. 1). Each quadrant of the mechanospace represents 
different combinations of flexural stiffness and torsional 
stiffness. ranging from the upper left, where beams twist 
easily but do not bend, to the lower right, where beams bend 
easily but do not twist. The dashed line in Figure | repre- 
sents a circular beam with a moderate Poisson’s ratio of 
0.25. The solid lines in Figure ] represent the two orienta- 
tions of a beam with a 4:1] elliptical cross-sectional shape as 
Poisson’s ratio is maximized (v = 0.5) or minimized (v = 
0). An elliptical beam of this form was chosen because such 
an ellipse has a large //J ratio and likely represents an 
extreme shape for biological structures. Thus, the area be- 
tween the two solid lines represents the range of values for 
beams with moderate to extreme material values, with 
shapes that are broadly representative of biological beams 
(Wainwright, 1988). Importantly, the solid lines do not 
represent absolute theoretical limits, but rather identify the 
expected extremes for ideal beams composed of a single. 


isotropic material. Greater shape modification will slightly 
alter these limits. 


Distribution of biological beams 


Although biological beams are commonly modeled using 
beam theory (e.g., Koehl, 1977b: Carrier, 1983; Vogel, 
1992; Baumiller, 1993: Ennos, 1993; Vogel, 1995; Niklas, 
1998; Etnier and Vogel. 2000; Etnier. 2001), they are rarely, 
if ever, made up of a single. isotropic material. More typi- 
cally, they are made up of multiple materials, whose distri- 
bution varies both across the cross section of the beam and 
along its length. Yet, in practice, flexural stiffness and 
torsional stiffness are measured experimentally using basic 
engineering formulas for beams. In general, a known load is 
applied to a beam, causing it to deform, either by bending or 
by twisting. The exact equation used is dependent on how 
the beam is loaded. For example, for end-loaded cantilever 
beams (one end fixed and the other end free to deform), 
flexural stiffness is calculated as: 


Di 


Sey 


(3) 
where F is the applied force. L is the length of the specimen. 
and y is the deflection of the free end of the specimen (Gere 
and Timoshenko, 1984). Similarly, torsional stiffness is 
calculated using the following formula: 


Fd 


GJ = oF 


(4) 
where F 1s the force applied at a moment arm d, L is the 
length of the specimen, and @ is the rotation of the free end 
in radians (Gere and Timoshenko, 1984). 

Data were compiled from the literature and from the 
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Figure 1. The mechanospace defined by values of flexural stiffness and torsional stiffness. The lines 
represent the distribution of ideal beams based upon an assumed cross-sectional shape and Poisson’s ratio. The 
symbols represent the nine groups (Table 2), with each individual point representing an average value for a 
species (appendix), The upper line represents an elliptical beam (4:t major to minor axis) with J// calculated 
about the short axis and a Poisson’s ratio of r = 0.50. The lower line represents that same beam with /// 
calculated about the long axis and a Poisson’s ratio of v = 0. The dashed line represents a circular beam with 


aePorsson S ratio of v = 0.25. 


author's research to obtain average flexural stiffness and 
torsional stiffness measures for 25 species of plants and 
animals (see appendix). These data were obtained following 
protocols similar to those discussed above. In other cases 
(n = 32), published data consisted only of measures of E 
and G. In these cases, the researchers experimentally deter- 
mined £7 and GJ, but then factored out / and J based on the 
cross-sectional size and shape of the structures. But note 
that the experimental data were still based on the overall 
mechanical behavior of the beam (i.e., deformation due to a 
load), rather than direct measures of material properties. In 
cases where published data consisted only of measures of E 
and G, both size and shape were estimated to determine 
values for 7 and J (appendix). Flexural stiffness and tor- 
sional stiffness were then calculated, based on these esti- 
mates. The estimates of beam diameter were deliberately 
conservative, while cross-sectional shape was always as- 
sumed to be circular. For example, for 14 of the 16 vines 
and 5 of the 11 tree branches included in the mechanospace, 
the beam was assumed to have a circular cross section with 
a diameter of 0.02 m. The reported diameters for these 
species were 0.02-0.05 m (Putz and Holbrook. 1991). Sim- 


ilar assumptions were made for the gorgonian corals, with 
an estimated diameter of 0.002 m (Jeyasuria and Lewis, 
1987). The assumption of a relatively small diameter will 
affect the overall magnitudes of flexural stiffness and tor- 
sional stiffness, but will not affect the ratio of the two. The 
assumption of a circular cross-sectional] shape will result in 
lower estimated values of //J than would be seen in noncir- 
cular beams (Table 1). 

The species (total n = 57) were divided into nine groups 
based upon broad morphological similarities (Table 2). 
Structures varied greatly in size, with diameters ranging 
from 0.001 m (red maple petioles) to 0.05 m (small tree 
branches). Average values for log GJ versus log El for all 
species (n = 57) were plotted in the mechanospace (Fig. 1). 
Note that each point represents an average value for a 
species, and thus does not reflect individual variation in size 
and material properties. For asymmetric beams, such as the 
daffodils, crinoid arms, and crustacean antennae, flexural 
stiffness is reported as the average of the different orienta- 
tions. Coefficients of variation for these organisms varied 
between 35% and 196% for flexural stiffness, and from 36% 
to 110% for torsional stiffness (Etnier and Vogel, 2000; 
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Table 2 


ed in this study 


Basic group characteristics for beams araly 


Groups ic characteristics 
‘osied beam. Antennae used for sensory 
ination and in some cases, for 

eressive interactions. Not loaded 
gravitationally. Aquatic animal. 

Multi-jointed beam. Arms used for passive 
filter feeding. Not loaded gravitationally. 
Aquatic animal. 

Multi-jointed beam. Stem must maintain 
upright position for photosynthesis, so self- 
supporting. Terrestrial plant. 

Continuons beam. Stem must maintain upright 
position for photosynthesis, so self- 
supporting. Terrestrial plani. 

Continuous beam. Petiole must support leaf 
against gravity and withstand wind, so self- 
supporting. Terrestrial plant. 


Crustacean antennae NI 


Crinoid arms 


Horsetails 


Herbaceous stems 


Leaf petioles 


Tree roots Continuous beam. Root is not self-supporting. 
Loaded in tension. Terrestrial plant. 
Vines Continuous beam. Vine is not self-supporting. 


Loaded in tension. Terrestrial plant. 
Continuous beam. Support individnal polyps 
against water flow. Not loaded 
gravitationally. Aquatic animat. 
Continuous beam. Tree must maintain upright 
position for photosynthesis, so self- 
supporting. Terrestrial plant. 


Gorgonian corals 


Tree branches 


The biological beams investigated were divided into nine groups, based 
on broad differences in morphology and function. 


Etnier, 2001), reflecting the individual variation noted 
above. 

The twist-to-bend ratio for each species was calculated 
and compared to the predicted values as an additional de- 
scriptor of flexibility in biological beams (Fig. 2). Standard 
deviations for the twist-to-bend ratios, when available, are 
reported in the appendix. 


Size-normalized mechanospace 


Size greatly influences the stiffness of a biological 
beam, both in twisting and in bending. The structural 
variables, / and J, are both proportional to radius to the 
fourth power; thus, small increases in size will greatly 
increase beam stiffness, regardless of cross-sectional 
shape or material composition. Values of E/ and GJ were 
normalized for size by dividing by radius to the fourth 
power to determine if size alone was driving the observed 
patterns. These values were then mapped into a size- 
normalized mechanospace (Fig. 3). This normalization 
accounts for size alone, rather than equating to calcula- 
tions of E or G for a given beam. 


Results 


The lines shown in Figure | are based upon the assump- 
tions of engineering beam theory and represent values for 
ideal beams. The distribution of biological beams closely 
matched that of the ideal beams, with 93% (53 of 57) of the 
points falling within the bounded region (Fig. 1). The 
boundaries reflect possible extremes for biological beams, 
based on assumed material values and chosen cross-sec- 
tional shapes. Flexural stiffness and torsional stiffness 
changed concurrently in the biological beams. Thus, beams 
that bent easily also twisted easily. while those that were 
hard to bend were also hard to twist. Overall, flexural 
stiffness and torsional stiffness each varied over 9 orders of 
magnitude. Species within the defined groups occupied sim- 
ilar regions of the mechanospace. implying that the flexural 
stiffness and torsional stiffness of group members were of 
similar magnitude. 

Due to the finite nature of the data set, the unoccupied 
regions of the mechanospace suggest that other combina- 
tions of mechanical values are rare without proving that 
they are impossible. Biological beams may exist within 
these unoccupied regions, given sufficient variation in ma- 
terial or structural design. The four samples falling outside 
of the predicted boundaries for ideal circular and elliptical 
beams included two herbaceous stems (daffodils Narcissus 
pseudonarcissus and sedges Carex acutiformis), a tree 
branch (Dendropanex arboreus), and a tropical vine (Marc- 
gravia rectiflora). These systems were all characterized by 
high flexural stiffness relative to torsional stiffness—that is, 
they twisted more easily than they bent (Fig. 2). 

Twist-to-bend ratios varied dramatically among different 
species (Fig. 2), with both the maximum (52.9) and the 
minimum (1.4) occurring in tropical vines. The average 
twist-to-bend ratio for all groups combined was 7.2, which 
falls between the predicted values for ideal beams that are 
circular or elliptical in cross-sectional shape. 

When the data were normalized for size (Fig. 3), the 
observed distributton changed slightly. In particular, the 
relative position of the antennae, horsetail rushes, and gor- 
gonian corals shifted up and to the right relative to the other 
beams. This result suggests that these structures are all 
relatively stiff for their size. While the position of groups 
changed slightly, structures that were stiffer in bending were 
also stiffer in torsion, regardless of their overall size. 


Discussion 


The close resemblance in the distribution of ideal and 
biological beams suggests that both are limited by the 
interdependence between the material properties, E and G, 
and between the geometric properties, / and J. Yet, keep in 
mind the assumptions of beam theory, which are violated 
almost universally by biological beams. The theoretical 
relationships apply only to beams consisting of a single, 
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Figure 2, Twist-to-bend ratios for the nine groups. The horizontal lines indicate the predicted values for 
ideal beams with circular (E//GJ = 1.00, v = 0) and elliptical (E//GJ = 12.75, v = 0.50) cross sections. The 
beams are arranged by size within each group, with the smallest diameter beam to the left. Standard deviations 
for the twist-to-bend ratios, when available, are reported in the appendix. 


isotropic material (Roark, 1943). Biological materials are 
more typically anisotropic, and values for Poisson's ratio 
can vary greatly from the theoretical expectations (Vincent. 
1990). More importantly. biological beams are almost al- 
ways composite structures built of multiple materials that 
differ greatly in their mechanical properties. Thus, it 1s not 
obvious that E/ and GJ would be highly interdependent in 
biological beams, particularly in the beams that deviate 
most from the theoretical assumptions. such as the jointed 
crinoid arms or crustacean antennae. 

The four samples falling outside of the predicted bound- 
aries are all characterized by high flexural stiffness relative 
to torsional stiffness—that is, they will twist more easily 
than bend (Fig. 2). The two herbaceous stems have flowers 
or seed heads that extend perpendicular to the long axis of 
the stem, potentially causing the stem to bend or twist when 
the wind blows. Rather than resisting this load with a high 
torsional stiffness, daffodils (Narcissus pseudonarcissus) 
and sedges (Carex acutiformis) have a low torsional stiff- 
ness, which allows them to twist in the wind. Daffodils 
reduce flow-induced drag with this action (Etnier and Vo- 


gel, 2000), and a similar function has been suggested for 
sedges (Ennos, 1993). The functional relevance of the high 
ratios of flexural stiffness to torsional stiffness in the tree 
trunk and vine has not been explored. 

Conservative estimates of size and shape were made for 
many of the tree branches, vines, and gorgonian corals 
(appendix). Beam size affects E7 and GJ equally; thus, 
errors in the size estimates will not affect the position of 
data points relative to the predicted boundaries (Fig. 4). 
Rather, an increase or decrease in diameter will cause these 
points to shift parallel to the boundaries. In contrast, shape 
estimates will differentially affect E7 and GJ. causing data 
points to shift relative to the predicted boundaries (Fig. 4). 
Although branches and vines are fairly circular in cross 
section, this assumption may not be valid for the gorgonian 
corals, whose cross-sectional shape can be circular, ellipti- 
cal, or even triangular (Jeyasuria and Lewis. 1987). The 
assumption of a circular cross section potentially underes- 
timates the range of values seen in the gorgonian corals. For 
example, if the cross-sectional shape of the gorgonians is 
assumed to be a 4:1 ellipse rather than circular, then all of 
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Figure 3. 


these beams will fall outside of the predicted boundaries. 
Future studies that include the cross-sectional shape of the 
gorgonian corals may identify group members that are of 
particular mechanical and functional interest within a bio- 
logical context. 

Flexural stiffness and torsional stiffness are both highly 
dependent on the size of the structure in question. While the 
pattern seen in the stiffness mechanospace is driven in part 
by size increases alone, the size-normalized mechanospace 
(Fig. 3) suggests a robust pattern. Even with differences in 
size removed from the comparison, there is still a strong 
relationship between how a beam bends and how it twists. 
The size-normalized mechanospace suggests some interest- 
ing comparisons between the different groups. For example, 
the antennae have a size-normalized flexural stiffness sim- 
ilar to that of tree branc} while exhibiting a higher 
size-normalized torsional stii . Again, the biological 
implications of such variation have not been sufficiently 
explored in the literature. 

The size, shape, and function of a beam may change 
during development (Wainwright and Dillon, 1969; Carrier, 
1983; Katz and Gosline, 1992; Gallenmiiller et al., 2001). 
Such changes potentially influence both flexural stiffness 
and torsional stiffness, thus affecting how the beam re- 
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Flexural and torsional stiffness normalized by size for all groups (Table 2). The data points were 
normalized by dividing by radius to the fourth power. 


sponds to a given load. The mechanospace introduced in 
this study is well suited for investigations relating the mor- 
phology, mechanics, and functional demands of biological 
beams during ontogeny. A systematic study of the relation- 
ship between EJ and GJ over a developmental series can 
reveal clues to the changing functional demands on a system 
during growth (Gallenmiiller et al., 2001). Concurrent mor- 
phological studies could determine whether mechanical 
changes are due to variation in material or structural prop- 
erties, potentially identifying the mechanisms that different 
systems use to modulate mechanical properties. 

Species within a group occupy a similar region of the 
mechanospace, signifying that the flexural stiffness and 
torsional stiffness of group members are of similar magni- 
tude. The same pattern is true for the size-normalized mech- 
anospace. While group members tend to occupy similar 
regions of the mechanospace, there is still variation within 
each group. Interspecific variation in shape and material 
may permit different combinations of flexural and torsional 
stiffness that may subtly influence how each species func- 
tions in its environment. 

Despite clear differences in design, no obvious mechan- 
ical differences distinguish jointed beams from continuous 
beams. In continuous biological beams, the observed me- 
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Figure 4. Changes in the assumptions of size and shape of beams will affect their distribution in the 
mechanospace. The small filled circles are the gorgonian values graphed in Figure 1 and were based on an 
assumed radius of 0.002 m. The lines represent the predicted values for ideal beams discussed previously. 
Changes in radius will affect the magnitude of 7 and J equally; thus, an increase or decrease in radius will cause 
these points to shift parallel to the predictions. For example, increasing the radius from r = 0.002 m to r = 0.004 
m shifts the distribution up and to the right (large circles). In contrast, shape estimates differentially affect Z and 
J. If the cross-sectional shape is assumed to be a 4:1 ellipse (fma, = 0.004 m, r,,,, = 0.001 m) rather than a circle, 
and 7 is calculated about the short axis, the data points move outside of the upper boundary (stars). Note that the 
cross-sectional area of the elliptical beams ts equivalent to that of the small, circular beams; thus, these changes 


do not reflect an increase in the total amount of material, but rather a change in its distribution. 


chanical properties may be attributed to materials that 
extend the entire length of the beam, but this explanation 
is inadequate for the jointed beams, which have no materi- 
als extending uninterruptedly along their entire length. 
The similarity in values for flexural stiffness and torsional 
stiffness for both jointed and continuous systems sug- 
gests that these beams may represent alternative designs to 
meet the functional need for flexibility in biological struc- 
tures. Two of the jointed systems, horsetails and antennae, 
are relatively stiff for their size, suggesting that the presence 
of joints does not necessarily equate with increased flexi- 
bility. 

Neither the ideal beams nor the biological beams are 
distributed uniformly throughout the mechanospace (Fig. 
1). Unoccupied regions of the mechanospace correspond to 
beams that bend but do not twis! and beams that twist bul do 
not bend. The distribution within the mechanospace may be 


determined by inherent principles governing the relation- 
ships belween E and /, as well as G and J. Conversely, 
emply areas within the distribution may not be an indication 
of physical impossibility, but of evolutionary history. There 
may be an absence of environmental patterns of change 
causing natural selection for particular combinations of me- 
chanical properties (Raup and Stanley, 1971). Alternatively, 
once an evolutionary pathway has been initiated, phyloge- 
netic canalization may limit future options for change 
(Lauder, 1982). Finally. empty spaces within the mechano- 
space may reflect temporal, rather than physical, limitations 
to those areas (Raup, 1966). The empty spaces will even- 
tually be occupied, given enough time. Although the iden- 
tification of boundaries within this mechanospace may not 
reveal their ultimate source, the boundaries do identify 
factors that may influence the observed pattern (Lessa and 
Patton, 1989), 
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Despite the structural diversity of the samples used in 
this study, they are merely a subset of the biological 
possibilities. A notable limitation is the absence of fiber- 
wrapped beams in the mechan:»pace. Internally pressur- 
ized, hydrostatic skeletons u pically wrapped with 
reinforcing fibers (Wainwright ez al, 1978). The fibers 
may be arranged ortho-onally with the fibers parallel and 
perpendicular to the long axis of the structure, or they 
may be in a helical array with fibers running in right- and 
left-handed helices around the long axis. Orthogonal 
arrays offer little resistance to twisting, while being rel- 
atively stiff in bending (Wainwright et al., 1978). Thus, 
these beams would potentially fall into the upper left- 
hand corner of the mechanospace. In contrast, helical 
fiber arrays allow pressurized beams to bend smoothly 
without kinking, while resisting torsional deforma- 
tions (Wainwright et al., 1978), potentially positioning 
these beams in the lower right of the mechanospace. 
Fibrous support systems may decouple the relation- 
ship between E/ and GJ, permitting novel combinations 
of mechanical properties. Thus, their inclusion in the 
mechanospace may greatly expand the observed distribu- 
tion. 

The mechanospace presented here is a useful approach 
for investigating patterns of flexibility in biological beams. 
Importantly, the mechanospace does not imply that flexi- 
bility has critical functional relevance in each system. 
Rather, it should be used as a broadly comparative tool to 
highlight systems in which flexibility may be biologically 
important. Biological beams that do not follow the basic 
pattern seen in the mechanospace may be of particular 
interest to both biologists and engineers, either due to ma- 
terial or structural innovation. 
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Appendix 


> along with their source, diameter, flexural stiffness (ED, torsional stiffness (GJ), and the ratio EI/GJ. 
cut their size and shape are marked with an asterisk. The standard deviations for E//GJ are given in 


Species (n = 57) used in the mechanospe: 
Species for which assumptions were ma: 
parenthesis, when available 


Diameter Ei GJ 
Group Species (m) (Nm’) (Nm‘) ENGJ Source 

Crustacean antennae n = 2 Procambarus sp. Crayfish 0.0004 2 30E—05 ~ G00E—06 4.5 (3.7) Etnier, 2001 
Panulirus argus Lobster 0.002 5.83E-03 2.04E-03 1.8 (1.3) Etnier, 2001 

Crinoid arms n = 2 Comactinia echinoptera 0.002 S00E [05 5" s0E— U5 4.3 (3.8) Etnier, 200) 
Florometra serratissima 0.002 3.92E—04  5.10E-—05 6.6 (2.7) Etnier, 2001 

Horsetails n = | Equisetum hyemale 0.005 PS7E S0 -3.83E—03 4.1 (1.8) Etnier, 1999 

Herbaceous stems n = 6 Cucumis sativus Cucumber 0.004 8.73E—03 oor US 5.4 (1-2) Vogel, 1992 
Helianthus annuus Sunflower 0.005 7.30E—05 5.30E-—05 1.4 (0.4) Vogel, 1992 
Lycopersicon esculenium Tomato 0.005 20E SO TII KGE = 03 ZI Vogel, 1992 
Carex acutiformis sedge 0.005 7.70E-03 2.34E-04 36.0(11.3) Ennos, 1993 
Tulips 0.006 20O E= 025. 2 10E 3 6.3.4.2) Etnier and Vogel, 2000 
Narcissus pseudonarcissus Daffodil 0.007 L.19E-02 8.90E-—04 13.3 (1.0) Etnier and Vogel, 2000 

Leaf petioles n = 5 Acer rubrum Red maple 0.001 1.94E-04° 0E =OS oie) Vogel, 1992 
Liquidambar styraciflua Sweet gum 0.002 9.84E-04 1.99E—-04 SAS) Vogel, 1992 
Phaseolus vulgaris Green bean 0.002 6.77E-04 1.43E-04 rke ah Vogel, 1992 
Populus alba White popular 0.003 T I3E=03 1. 50E=05 4.95 Vogel, 1992 
Populus tremuloides 0.003 KOQE =05 — 1.00E 506 4 Niklas, 199} 

Tree roots n = | Pinus taeda Loblolly pine 0.05 I 72EF02 7 90E7 01 2.3 (0.8) Vogel, 1995 

Vines n = 16 Croton pullei juvenile 0.01 3.30E+00 3.04E—-0]1 10.9 Gallenmiiller er al., 2001 
Cissampelos pareira 0.023 3.63E+00 4.90E—O1 7.4 Putz and Holbrook, 1991 
Cissus sicyoides O02" RISESOGO™ 2)30E—C1 8.8 Putz and Holbrook, 1991 
Forsteronia portoricesis 0.02* 1.67E+00 8.30E—0l1 2.0 Putz and Holbrook, 1991 
Heteropteris laurifolia O024 107E+01 4.60E+00 2 Putz and Holbrook, 1991 
Hippocratea volubilis 0.02* 4.38E+00 6.80E—-01 6.4 Putz and Holbrook, 1991 
Ipomoea repanda 0027 JAIE FOOT 3.1CE—01 11.0 Pntz and Holbrook, 1991 
Marcgravia rectiflora 02. 1.32E FO less 0E—01 JAS Putz and Holbrook, 1991 
Mikania fragilis 0.02* 1.48E+00 2.00E-01 7.4 Putz and Holbrook, 1991 
Paullinia pinnata 0.02* 471E+00" 6.30E—01 7.5 Putz and Holbrook, 199] 
Rourea surinamensis 0.02* 648E+00 1.23E+00 5.3 Putz and Holbrook, 1991 
Schlegelia brachyantha 0.02* 4.34E+00 9.60E—0]1 4.5 Putz and Holbrook, 1991 
Securidaca virgata O02? 1.59E+60 3.60E—01 4.4 Putz and Holbrook, 199} 
Croton pullei mature 0.02* 4.87E+00 » 3.616700 1.4 Gallenmiiller er al., 200} 
Vitus rotundifolia Grape 0.046 3.08E#02  LI4EFO? 25) (03) Vogel, 1995 
Wisteria sinensis Wisteria 0.051 283E t0? E6 T0E FOI 4.5 (2.1) Vogel, 1995 

Gorgonian corals n = 13 Ellisella barbadensis 0.004* KiGE=—01_ "2° 16E—02 53 Jeyasuria and Lewis, 1987 
Eunicea calyculata 0.00-4* Oj E03 AI OES) Sl Jeyasuria and Lewis, 1987 
Eunicea clavigera 0.004* TWE 03 2-26E—03 35 Jeyasuria and Lewis, 1987 
Gorgonia ventalina 0.004* 160E- 0 = IO0S5E=02 1.6 Jeyasuria and Lewis, 1987 
Leptogorgia virgulata 0.004* 1I7E = 0? IBE 02 3/3 Jeyasuria and Lewis, 1987 
Lophorgogia cardinalis 0.004* 635E -02 — 1 46E—02 4,3 Jeyasuria and Lewis, 1987 
Muriceopsis flavida 0.004* Ooch ao 243. IE 03 2 Jeyasunia and Lewis, 1987 
Plexaura flexuosa 0.004* 1256-702 . 2 Ole U3 62 Jeyasuria and Lewis, 1987 
Plexaurella grisea 0.004* 5.66E—0: TIME- 03 5.9 Jeyasnria and Lewis, 1987 
Pseudoplexaura crucis 0.004* T3J1E -0 T 2OIE -03 3.6 Jeyasuria and Lewis, 1987 
Pseudopterogorgia bipinnata 0.004* 4172-020" ZOPET 03 4.6 Jeyasuria and Lewis, 1987 
Pterogorgia citrina 0.004* LISE = 0 ys UIE-05 6.5 Jeyasuria and Lewis, 1987 
Swiftia exserta 0.004* 8. 82E—03 1.76E—03 5.0 Jeyasuria and Lewis, 1987 

Tree branches n = 1] Brunellia comocladifolia 0.024 2.08E+01  1.65E+409 12.6 Putz and Holbrook, 1991 
Dendropanaa arboreus 0.02* 9.42E+01 2.81E+00 235 Putz and Holbrook, 1991 
Guarea trichilioides 002% LOJE+01 1.95E480 9.9 Putz and Holbrook, 1991 
Inga vera 0.02* 3.74E+01 5.47E+00 6.8 Putz and Holbrook, 1991 
Ocotea floribunda 0.02* IT3JIEFOI 4.23E+09 Bel Putz and Holbrook, 1991 
Juniperus virginiana E. red cedar 0.04 7.44E+02 1.68E+02 4.4 (0.7) Vogel, 1995 
Liriodendron tulipifera Tulip poplar 0.045 LIGEt03 1 3002 8.9 (1.4) Vogel, 1995 
Acer rubrum Red maple 0.047 2.072403 250E+02 8.3 (0.9) Vogel, 1995 
Phyllostachys sp. Bamboo 0.05 8.03E+03 9.50E+02 8.6 (1.3) Vogel, 1995 
Pinus taeda Loblolly pine 0.051 206E t03 3.19E+02 6.1 (0.9) Vogel, 1995 
Liquidambar styraciflua Sweet gum 0.054 233E F03 2.60E702 8.9 (1.6) Vogel, 1995 


